DOE/ER/40561-296-INT96-00-153 

UW/PT-96-31 

nucl-th/9610052 



More Effective Field Theory 
for Nonrelativistic Scattering 



David B. Kaplan 
Institute for Nuclear Theory, University of Washington 
Box 351550, Seattle WA 98195-1550 
dbkaplanOphys . Washington . edu 



An effective field theory treatment of nucleon-nucleon scattering at low energy shows 
much promise and could prove a useful tool in the study of nuclear matter at both ordinary 
and extreme densities. The analysis is complicated by the existence a large length scale 
- the scattering length — which arises due to couplings in the short distance theory 
being near critical values. I show how this can be dealt with by introducing an explicit 
s-channel state in the effective field theory. The procedure is worked out analytically in a 
toy example. I then demonstrate that a simple effective field theory excellently reproduces 
the 1 5'o np phase shift up to the pion production threshold. 
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1. Introduction 

Effective field theory is a powerful tool in particle physics for describing a scatter- 
ing amplitude in a physical system with several length scales (for recent reviews, see [II). 
It has been used with great success to describe both relativistic systems (e.g., the Stan- 
dard Model of weak interactions) and nonrelativistic ones (e.g., positronium 0, as well as 
bound states containing heavy quarks HU). Such applications involve a systematic ex- 
pansion in external momenta that allows one to compute low energy physical observables 
consistently to a desired accuracy. Recently a number of papers have applied effective 
field theory to nucleon-nucleon interactions, with an eye toward better understanding the 
properties of nuclei and nuclear matter |5|-|8|]. In this case one is dealing with nonrela- 
tivistic scattering and a short-range interaction, but the power expansion of the effective 
theory is complicated by the presence of an independent length scale in the form of large 
scattering lengths. For example, np scattering in the 1 5'o channel has a scattering length 
a = — 23fm = — l/(8MeV). In this letter I show how to construct an effective theory for 
such a system. I demonstrate the method in an analytically soluble model, and then apply 
it to the problem of 1 5'o np scattering. I find that the method, which includes in the 
effective theory a fundamental field exchanged in the s-channel, describes np scattering far 
more successfully than previous effective field theory calculations (see fig. 5 below). 

To motivate what follows, I briefly review the justification for effective field theory. 
Experiments at momenta p <C A are insensitive to the effects of cuts and poles in the 
amplitude associated with the heavy states at p > A. Therefore the amplitude can be well 
approximated by a function which correctly reproduces the nonanalytic features associated 
with light states at scales p C A, plus a smoothly varying function which can be expanded 
in powers of p/ A to account for the effects of the distant heavy states. This can be achieved 
by introducing an effective Lagrangian £ e ff that includes the light states as propagating 
degrees of freedom, along with nonrenormalizable local interactions in a power expansion in 
derivatives and fields suppressed by the appropriate powers of A. The masses and couplings 
of the light states are adjusted to put the low lying singularities of the amplitude in the 
right places. 

If the correct short distance theory is known, this matching of low lying singularities 
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can sometimes be done by explicitly by "integrating out" the heavy degrees of freedom. In 
this case, the effective theory is a tool that greatly simplifies calculations by eliminating 
degrees of freedom that decouple from low energy phenomenology. 

On the other hand, if the correct short distance theory is not known, or the matching 
cannot be performed, then the couplings in the low energy theory are determined from 
data. An example of this "bottom-up" effective theory is the chiral Lagrangians for pion 
physics. Such theories typically are able to describe low energy physics in terms of few 
parameters, and are very predictive. 

This letter is organized as follows: I first perform a top-down calculation for a soluble 
system with short distance interactions (the three dimensional 5-shell potential), perform- 
ing the nonperturbative matching onto the effective theory analytically, showing how to 
avoid problems as the interaction is tuned such that the scattering length diverges. I then 
generalize the matching procedure to a system with both short and long distance interac- 
tions. I then apply the method with great success to np scattering in the 1 Sq channel. The 
latter is a bottom-up calculation, where the coefficients of the effective theory are chosen 
to best fit the data. For simplicity I restrict my discussion to s-wave scattering throughout 
the paper. 

2. A toy model: the 5-shell potential 

It is instructive to examine how the effective field theory works for an analytically 
soluble short distance interaction. The example I use is of two "nucleons" of mass M 
interacting in three dimensions via the potential V(r) given by 

V(r) = -g±6(r-l/A) , (2.1) 

where r = |ri — r 2 | is the internucleon separation. The parameter A has dimensions of 
mass (I set c = h = 1) and sets the size of the interaction region. The coupling g is 
dimensionless, and for s-wave scattering, there is a single boundstate for g > 1 and no 
boundstate for g < 1. 

In the center of mass frame of the two nucleons, the s-wave Feynman scattering 
amplitude A = —iT is given in terms of £ = p/A = \[ME / 'A and the spherical Bessel 
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functions jo and no 



A - ^ ( ) (22 ) 

This amplitude has a simple pole at 

1 - e~ 2v 1 

£ = iv, where — = - ; (2.3) 

2r] g 

and |£o| — *■ as g — > 1. This singular behavior is evident from the scattering length a 
computed from (|27|). One finds 

9 ^ 1 (2.4) 



1-gJ A ' 

which diverges as g — > 1. Evidently a naive expansion of the amplitude about p < A 
(small £) will fail for |^| > — g)/g\. It is the dimensionless factor of g/(l — g) that makes 
simple dimensional analysis for this system unreliable for g ~ 1. 

Nevertheless, all other poles in the amplitude ( |2.2| ) occur at values |£| > 2%. Following 
the precepts of effective field theory, the reasonable approach is to describe this system 
as one of nucleons N, a fundamental field (p which has baryon number B = 2 and mass 
2M + A (henceforth called "the dibaryon" ) , and local interactions expanded in powers of 
V/A. The field is included to reproduce the low-lying pole in the amplitude; a radius of 
convergence of \p\ > A is expected for the residual interactions, since all other singularities 
in the amplitude are at that scale. The effective Lagrangian is given by 

£ eff =iV t Ud t - V 2 /2M) N + atf (id t - V 2 /4M - A) 

(2.5) 

- (ypNN + h.c.) - \C(N^N) 2 + ... 

where the ellipses refer to higher dimension operators suppressed by powers of A. The field 
(ft can be rescaled so that \a\ = 1, but the sign of a will have to be determined. In the 
center of mass frame, the scattering amplitude can be calculated by summing the ladder 
diagrams (solving the Lippman-Schwinger equation) with the kernel 

7eff(p)=Cf+ ^TA' E = p 2 /M, (2.6) 
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Fig. 1. The potential from the effective theory ( \2.3j ), and the amplitude j \2.8j) ob- 
tained from summing the ladder (bubble) diagrams. 

as shown graphically in fig. 1. The bubble diagrams in the geometric sum are divergent, 
but renormalization of C, A and y renders the theory finite. Of course, the values for these 
renormalized parameters are scheme dependent and not physically meaningful. What is 
desired is a scheme that facilitates the power counting and effective field theory expansion 



in p/A. Dimensional regularization and minimal subtraction (MS) is such a scheme, as 
the renormalization scale \i only appears in logarithms; I use this scheme throughout the 
paper. The basic calculation one needs for the present model is the bubble integral M 



d 3 q 1 MV-ME - it -iM\p\ 



(2tt) 3 (E - q 2 /M + ie) An An 

The resultant amplitude is 



(2-7) 



„ U - YM - C(aE-A) + y 2 

l + iMV(p)p/Air aE - A + i(Mp/An)[C(aE - A) + y 2 } ' v ' 



In eqs. ( |2.2|) and ( |2.8| ) I have written the exact and effective amplitudes in the form 
N(p)/D(p) where iV and D are nonsingular functions of p. By correctly choosing the 
parameters of the effective theory, it is possible to match iV and D in the two theories 
up to 0(p 4: ) corrections and a common multiplicative constant which factors out of the 
amplitude. This involves solving the following equations to order 0(p 2 ): 

-[C(^-A)+^]=*i^fc«)]S, 

(aE-A) = z(l+gOo(Z)Mt)) , 
where z is the irrelevant constant. 

One finds the following solution to the matching equations 



(2.10) 



A = -a (^) = a^- , f = - (111) ^ 

V g J 2M oM' y \ g J M 2 



The sign of a is determined by the fact that y 2 must be positive: 



a = -Sign . (2.11) 

There are several interesting features about these values: 

(i.) The parameter A is proportional to the inverse scattering length 1/a, and goes to zero 
at the critical coupling g = 1. Since the momentum variation of the (p propagator 
depends on the combination p 2 /MA, for g ~ 1 the </> field is needed in the effective 
theory to reproduce the rapid variation of the scattering amplitude at low p. On 
the other hand, for \g/(l — g)\ ~ then the scattering length is short and the 

4> field looks heavy (MA ~ A 2 ); it can be integrated out without compromising 
the convergence of the momentum expansion, but then the effective theory is simply 
perturbative 0. 

(ii.) The parameter C is independent of g and is small, in the sense that CMp/ (4tt) ~ p/A, 
so that its effects only become nonperturbative for large momenta, p ~ A. 

(iii.) For a = —1, which is the case when \g + h\ > |, the <fi field in the effective theory 
( p . 5|) has the wrong-sign <9 t energy term. This is no cause for alarm, since the nonper- 
turbative amplitude one computes using the fully dressed 4> propagator and vertices is 
unitary and sensible. Loops with internal <p particles are not included in the effective 
theory. 

In fig. 2 I compare the exact phase shift from eq. ( |2.2|) with the effective field theory 
prediction ( |2.8| ) for the two values g = 1 ± .01. These are systems where the scattering 
length a in eq. (|2.4| ) is 100 times the fundamental length scale of the problem, 1/A. As one 
can see, the size of the factor g/{l — g) does not affect the range of validity of the effective 
field theory. In fig. 3 I do the same for very strong couplings g = ±10. As expected, in 
every case one compares, whether it be strong or weak coupling, attractive or repulsive, 
the effective theory agrees well with the exact result for momenta p < A. Evidently the 
effective theory ( |2.5| ) is very effective at accurately describing long distance physics. 

If one tries to do without the dibaryon state cf) in the effective theory, using instead 



1 For an enlightening discussion, see ref. [10]. 




Fig. 2. s-wave phase shift (degrees) for the 5 -shell potential, plotted as a function of 
p/A for g = 0.99 (with an almost-bound state) and g = 1.01 (with a weakly bound 
state). The solid lines are the exact results, while the dashed lines are computed 
from the effective field theory as described in the text, with M/A = 1.2. 




Fig. 3. s-wave phase shift (degrees) for the d-shell potential, for g = 10.0 (strongly 
attractive) and g = —10.0 (strongly repulsive). The solid lines are the exact results, 
while the dashed lines are computed from the effective field theory with M/A = 1.2. 



the effective Lagrangian 

C eS = jV f (id t - V 2 /2M) N - \C{N ] N) 2 - \C 2 (N ] VN) 2 + ... (2.12) 
then one runs into difficulties. The coupiings of this theory £ e ff can be readily obtained 



from C e ff in eq. (|2.5|) by integrating out the <fi field. This amounts to performing a 
Taylor expansion of interaction V e ^{p) in eq. ( |2.6| ) for small p. One finds C = 4na/M 
(as was found in ref. ||), and that higher derivative operators in £ e ff form an expansion 
in powers of p 2 /MA ~ p 2 a/A. This expansion is fine so long as there is no boundstate 
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near threshold. However, when a is large, A is small and £ e ff does not constitute a useful 
expansion. Such is the case for realistic s-wave NN scattering, for which the scattering 
lengths are much larger than the range of the short distance interaction. 

This observation exemplifies the general discussion of ref. [Q: that the radius of 
convergence of the derivative expansion an effective field theory without a dibaryon is set 
by the scale p ~ a/ A/a, which is poor for systems with large scattering length a. However 
it was argued in || that in fact the higher derivative operators were all correlated, and that 
the expansion of 1/A would have a radius of convergence p ~ 1/A. The method described 
here of including an explicit dibaryon field in the s-channel is a more straightforward and 
intelligible procedure containing the same physics. 



3. A toy model with both short- and long-range interactions 

Performing an effective field theory calculation for nucleon interactions is complicated 
by the existence of both short- and long-range interactions. The latter is due to one 
pion exchange (OPE), which must be included explicitly in the effective theory. One 
might worry that the effects of pion interactions could completely change the nature of the 
solution ( |2.10| ), particularly if pion exchange is strong enough to contribute significantly 
to the total scattering length. It is therefore instructive to complicate the 5-shell model 
of the previous section by including a Yukawa interaction with range 1 / ^> 1/ A. Using 
the techniques developed in ref. 0, it is possible to perform the matching analytically. 
In fact, I will show that even if the pion interaction is strong, its effects on the matching 
coefficients is computable and small — (9(^ L ). The reader uninterested in technical details 
can skip to §4. 

3.1. The amplitude of the full theory 

For the interaction in the full theory I take the potential V(r), where 

V = V s + Vir, V s = -^8(r -1/A) , V* = -^e~ m - r . (3.1) 

M r 

Following ref. [|], I define the OPE Hamiltonian = — V 2 /M + V n with dimensionless 
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solutions J p , Mb satisfying i 

(H ir -E)J p (r) = 
with expansions near the origin 



(H n -E)Af p x (r) 



4n 



Mp 



(3.2) 



J p (r) -> 1 - ^Lr + ^M + alM*-2 P * r2 + ^ 



+C(r 2 In Ar) 



r (3.3) 



The parameter A is due to the ambiguity of being able to add the solution J to J\f while 
still satisfying eq. ( |3.2| ); nothing physical will depend on it. 

An s-wave scattering solution with momentum p to the Schrodinger equation with 
the full potential V given in eq. (|3.1|) can be written as tp p (r) = aj p {r) + bj\fp(r) for 
r > 1/A, while ifj p (r) oc J p {r) for r < 1/A. It is straightforward to compute the ratio a/b 
by performing the matching at r = 1/A, making use of the Wronskiani: 



a _ M 
b 4n 



-P + 



A 



(3.4) 



=1/A 



The s-wave phase shift is then determined by a/b and the asymptotic properties of J 
and M . The latter can be expressed as 



J P (r) — ► { r qe Lpr jpr + c.c.) 



((e ipr /pr + c.c.) 



(3.5) 



The phase shift is then given by 



(ay + b() 
(arj + bCY 



(3.6) 



2 Comparing with analogous functions defined in ref. 
-(4iv/Mp)K E (r), where p = \[EM. 

3 W — Kf x — 7 — 7 — \f x 



I, Mr) = Je{t) and M^r) 



l 

p r 2 
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where 77, £ are the asymptotic coefficients of eq. ( |3.5| ). Following ref. Q, it is convenient 
to use the expressions 



77/77* = e 2 ^ , 

I/77* = 27 Xp (0) , (3 7) 



C7^ = Gf s (o,o) + ^ 



In the above expression, 5 n , Xp(0), Gg S (0, 0) are quantities computed from the OPE 
Hamiltonian H n : 5 n is the OPE phase shift, x p (0) is the OPE wavefunction at the ori- 
gin, and Gg S (0,0), the Green function (r| 1/(12 — H n ) \r') at r' = r = 0, dimensionally 
regulated and renormalized in the MS scheme, at the renormalization scale \x @. 

Combining eqs. (|3.4|) - (|3.7| ) I determine the complete scattering amplitude to be 



t^4. */4-7j- I 



\MA J 



1 + g&TfJp + 9^kJ p 2 (in $ + 2 7 - l) + G% s (0, 0) 



In this expression, J7" p and Af£ are evaluated at r = 1/A. 

Note that in the limit that the pion coupling vanishes (a^ — > 0) the terms in the above 
expression behave as 



At — ► , x P (0) — 1 , Gf s (0,0) — > - ^ , (3.9) 
while the J7" and M functions become spherical Bessel functions: 



Mr) ~^Mw) , AT A (r) — » n (pr) . (3.10) 

Ott— >0 a„-— >(J 

Upon substituting the above expressions in eq. Q3.8TJ one recovers the amplitude ( |2.2|) for 
the 5-shell potential alone, a T = 0. 



The calculation of the amplitude A never requires renormalization; I express the answer in 
terms of the renormalized propagator to more easily compare to the effective theory, which does 
require renormalization. 
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V eff= U. = X + > """" < + 




7C 





Fig. ^. TTie potential from the effective theory including the explicit dibary on field 
cp and single pion exchange. The amplitude A e ^ is obtained from summing the 
ladder (bubble) diagrams. Techniques for performing the sum using dimensional 
regularization are extensively discussed in ref. f^j. 

3.2. The amplitude in the effective theory 

The effective theory with which I expect to describe fow energy scattering consists of 
a contact interaction and both <p and 7r exchange. The potential for the effective theory 
is given graphically in fig. 4; it is the same as the effective interaction in fig. 1, with the 

addition of the single pion exchange graph: 

2 

v a 

V eS = V + V n , V Q (E) = C+ y , K(p,P , ) = - 7 T^T 2 • (3-11) 

aE — A (p — p')^ + 

Using the techniques of ref. || it is possible to write the resultant scattering amplitude in 

closed form as 

A V (E)( Xp (O)f 



l-U (F)Gf s (0,0) 

[C(aE-A) + y 2 ] ( Xp (0)f 
aE — A - Gf 5 (0, 0) [C(aE - A) + y 2 } 



(3.12) 



3. 3. Matching the effective and full amplitudes 

Matching entails choosing the free parameters {y 2 , A, C} such that the amplitude 
( |3.12|) best approximates the exact exact amplitude (|3lf). The analogue of the matching 



equations (|2.9|) in terms of £ = p/A are: 



4-7T 

■[C(orE-A) + y 2 ]=z m g[J p ] 2 



(aE - A) 



1 + gZJrtfi + 9^J P 2 (in ^ + 2 7 - 1 
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(3.13) 



where again J p and A/3 are evaluated at r = 1/A. These matching conditions are to be 
solved to 0(p 2 ). 

So what has been the effect of the pions? The effective theory is supposed to be 
an expansion in p/A; must one perform a simultaneous expansion in a n 7 To resolve this 
question note that one pion exchange must be treated nonperturbatively for a 7T M/m n > 1. 
In the real world, this value is nearly realized — for NN scattering in the 1 5'o channel, 
a n = 0.47m n /M. This suggests that one should consider a n ~ m^/M in any power 
counting scheme. Furthermore, for momenta < p <C A, the effective theory is an 
expansion in e = p/A > m^/A. Combined, these considerations imply 

p a^M 

A A A (3.14) 

a^Mm^ a 2 M 2 2 

A 2 A 2 



Thus when J v {r) is evaluated at r = 1/A for the matching condition (|3.13|) , I find from 
eq. O 

J p {l/A) = (1 + 0(e)) - (1 + 0(e)) . (3.15) 

An important lesson has been learned from this exercise: one can perform the matching 
in a perturbative expansion in e ~ m w /A, even when pion interactions are nonperturbative 
and a n M/m 7T ~ 1. Since a n enters the matching via the combinations in eq. ( |3.14j ), the 
leading order calculation in e corresponds to setting a n to zero, and one recovers the 
solution (p.lO| ) found when performing the matching without pions. At subleading order, 
the 0{e) contributions from matching are more important than including two derivative 
operators in £ e ff, which are 0(e 2 ) effects. One may be tempted to think that this is 
a manifestation of chiral symmetry, since a n oc m^/ f 2 and vanishes in the chiral limit. 
However, this is not the case, since we are assuming that a n M/m 7T ~ O(l), while that 
combination would vanish in the chiral limit. 

Finally note that the In terms in Af£ combine with the In term in eq. (|3.13[ ) to 

2 

give In j^. This logarithm accounts for running of the effective interactions at scales \x 
below the matching scale A. 
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4. NN scattering in the 1 Sq channel 



I now turn to a problem of phenomenological interest: np scattering. As discussed 
above, this is a system with a scattering length a = —23 fm much larger than the expected 
scale of short-range strong interactions (which are characterized by a scale shorter than 1 
fm). Therefore an effective field theory description can be expected to benefit by explicit 
inclusion of a dibaryon field <f> as discussed in the previous sections. The effective theory for 
np must also include explicit pion fields. Since the short-distance theory is not understood, 
I fit the couplings of the effective theory to data, rather than computing them via matching 
conditions. The results are expected to be similar to those in ref. ||, where the <fi field 
was not included explicitly, but where 1/A was expanded so as to be insensitive to the <p 
pole. One advantage of the present method is that unlike the treatment in ||, the present 
method should be valid even for small phase shift (where 1/A blows up). Furthermore, 
it is possible to account for the residual interactions in an expansion in p/A, where A 
characterized short distance physics. 

The potential summed for np scattering is given graphically in fig. 4 with 

V eS EE V + V n , 

(p - p ) + m i 

The pion axial coupling g& = 1-25 and decay constant = 132 MeV are known; thus the 
effective theory involves the same three unknown parameters as in the 5-shell example: 
The four-nucleon interaction C, the dibaryon mass shift A, and the dibaryon coupling y 
(as well as the sign a of the dt kinetic term). As before, the values for these quantities are 
not physical, depending on the renormalization scheme chosen. The scattering amplitude 
one derives, however, is independent of the renormalization scheme. 

The amplitude A e ff of the effective theory is the same as appearing in eq. ( |3.12| ), 
with Vq and a n given by eq. (|4.1|) . By numerically computing the quantities S n , x an d G, 
and subsequently fitting the resultant phase shift to the data, I find the following (scheme 
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(4.1) 



9A m l 



100 200 300 400 

P Km (MeV) 



Fig. 5. 1 Sq phase shift (degrees) versus center of mass momentum p cm for n — p 
scattering. (p cm = 400 MeV corresponds to = 340 MeV). The solid line is 
taken from the Nijmegen partial wave analysis j\Tl\J, while the dashed line is computed 
from the effective field theory in eq. (\3.12j) , which involves three free parameters. 



dependent) values for the unknown parameters at renormalization scale fx = m^: 



C = - 



1 



192 MeV 
1 



C 



c + 



9 2 A 



507 MeV 



A = 14.8 MeV 



1 



237 MeV 



(4.2) 



and a = — 1. Note that the contact interaction can be written as C — 47r/MA, where 
A = 749 MeV; thus a loop with a contact interaction brings with it a power of G ^ p = p/A. 
As this is independent of M, the power counting scheme is different in such a theory from 
that described in refs. all residual interactions in the derivative expansion of £ e ff can 

be treated perturbatively. 

As is shown in fig. 5, the phase shift 5 e s computed from this effective theory with 
its three free parameters is essentially indistinguishable from the results of the Nijmegen 
phase shift analysis for 1 5'o scattering, all the way up to center of mass momentum p cm = 
400 MeV, corresponding to Ti a b = 340 MeV. That the fit works so well is evidence that 
higher derivative are suppressed by powers of a large scale, while the relevant effects of 
multipion exchange, A exchange, etc. are well accounted for by the effective interactions 
included in £ e fj in eq. (|4.1| ). 
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5. Conclusions 

Expanding on techniques introduced in ref. ||, I have demonstrated that effective 
field theory and dimensional regularization can be successfully applied to the study of 
nonrelativistic systems with short-range interactions. I have shown that when the short- 
range interaction is near the critical value for which the scattering length diverges, then a 
low-lying s-channel state must be introduced as a physical degree of freedom in the effective 
theory. This procedure is in keeping with how effective field theory is applied to relativistic 
systems when a mass independent renormalization scheme is employed PJ. As shown by 
explicit example, the method works well for both attractive and repulsive interactions, 



despite the fact that problems can arise when using other forms of regularization ||12|| . 

The effective theory reproduces the nonperturbative features of interactions at the 
scale A through the s-channel resonance, plus local interactions which can be treated 
perturbatively in an expansion in p/A. Including additional long-range, nonperturbative 
interactions does not disturb this power counting. 

The method also works very well in reproducing np scattering data in the 1 So channel. 
It is hoped that the technique can be used for the spin triplet channels as well, but that 
has yet to be resolved. Eventually, the value of the approach outlined here is hoped to lie 
in the study of nuclear matter at both normal and extreme density. While the usefulness 
of effective field theory has yet to be determined for bulk matter, it has been shown 
here to reproduce the two particle scattering amplitude at well above the nucleon Fermi 
momentum in matter, so there is reason for optimism. 
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